Soyasapogenol-B as a Potential Multitarget Therapeutic Agent for Neurodegenerative Disorders: Molecular Docking and Dynamics Study

Neurodegenerative disorders involve various pathophysiological pathways, and finding a solution for these issues is still an uphill task for the scientific community. In the present study, a combination of molecular docking and dynamics approaches was applied to target different pathways leading to neurodegenerative disorders such as Alzheimer’s disease. Initially, abrineurin natural inducers were screened using physicochemical properties and toxicity assessment. Out of five screened compounds, a pentacyclic triterpenoid, i.e., Soyasapogenol B appeared to be the most promising after molecular docking and simulation analysis. Soyasapogenol B showed low TPSA (60.69), high absorption (82.6%), no Lipinski rule violation, and no toxicity. Docking interaction analysis revealed that Soyasapogenol B bound effectively to all of the targeted proteins (AChE, BuChE MAO-A, MAO-B, GSK3β, and NMDA), in contrast to other screened abrineurin natural inducers and inhibitors. Importantly, Soyasapogenol B bound to active site residues of the targeted proteins in a similar pattern to the native ligand inhibitor. Further, 100 ns molecular dynamics simulations analysis showed that Soyasapogenol B formed stable complexes against all of the targeted proteins. RMSD analysis showed that the Soyasapogenol B–protein complex exhibited average RMSD values of 1.94 Å, 2.11 Å, 5.07 Å, 2.56 Å, 3.83 Å and 4.07 Å. Furthermore, the RMSF analysis and secondary structure analysis also indicated the stability of the Soyasapogenol B–protein complexes.


Introduction
Alzheimer's disease (AD) is the most common form of dementia, characterized by a gradual deterioration in cognitive abilities [1]. Pathological features of the condition include amyloid-beta (Aβ) and tangled intracellular neurofibrils with aggregated tau protein (hyperphosphorylated), as well as synaptic and neuronal loss [2][3][4]. Currently, there are more than 55 million people affected by AD globally, with an increase of 10 million new cases every year, which is predicted to grow to 106.4 million by 2050 [5]. The rapidly increasing prevalence of AD is likely to put a major burden on families and the community [1]. olated Lipinski's rule, and only 18 compounds were found to have less than 90 Å TPSA (required to cross blood-brain barrier) and better (more than 70%) gastro-intestinal absorption (Supplementary Table S1). We further analyzed the toxicity potential of 30 BDNF inducers (those that did not violate Lipinski's rule), and found that only 13 compounds exhibited no toxicity for four parameters: mutagenic, tumorigenic, reproductive effect, and irritant. Among these 13 non-toxic compounds, we found that only five compounds were among those 18 compounds that were selected earlier on the basis of their intestinal absorption percentage and lower TPSA, suggesting their efficiency in crossing the blood-brain barrier. Soyasapogenol B, Huperzine A, Hydroxytyrosol, Alpinetin, and Calycosin showed good absorption efficiency without violating the Lipinski rule (Supplementary Table S1), with minimal or no toxicity (Supplementary Table S2).

Analysis of Interaction of Soy-B with AChE
In the study, interaction between Soy-B and AChE revealed that Soy-B interacted with the residues close to the catalytic active site of AChE, and hence showed a noncompetitive inhibition (Figure 1). The binding mode of Soy-B was further compared with native control ligand (Tacrine), and the protocol was standardized by re-docking the native ligand ( Figure 1B). Native ligand in the co-crystalized structure interacted with TRP84 and PHE330 amino acid residues of AChE via pi-pi stacking and pi-alkyl interactions, respectively. Similar results were observed when native ligand was redocked with AChE ( Figure 1E). Moreover, eleven amino acid residues (SER 81, GLY117, GLY118, GLU199, SER200, TYR 334, TRP432, ILE439, HIS440, GLY441, TYR442) of AChE were involved in van der Waals' interactions with redocked ligand. On the other hand, Soy-B interacted with AChE near to the native ligand binding site ( Figure 1C) through two strong hydrogen bonds via ASP285 and GLY335 amino acid residues of AChE. Stability of Soy-B and AChE binding was further enhanced by van der Waals' interactions with TYR70, ASP276, VAL277, TRP279, ASN280, SER286, ILE287, TYR334 and LEU358 amino acid residues ( Figure 1D).

Analysis of Interaction of Soy-B with AChE
In the study, interaction between Soy-B and AChE revealed that Soy-B interacted with the residues close to the catalytic active site of AChE, and hence showed a non-competitive inhibition (Figure 1). The binding mode of Soy-B was further compared with native control ligand (Tacrine), and the protocol was standardized by re-docking the native ligand ( Figure 1B). Native ligand in the co-crystalized structure interacted with TRP84 and PHE330 amino acid residues of AChE via pi-pi stacking and pi-alkyl interactions, respectively. Similar results were observed when native ligand was redocked with AChE ( Figure 1E). Moreover, eleven amino acid residues (SER 81, GLY117, GLY118, GLU199,  SER200, TYR 334, TRP432, ILE439, HIS440, GLY441, TYR442) of AChE were involved in van der Waals' interactions with redocked ligand. On the other hand, Soy-B interacted with AChE near to the native ligand binding site ( Figure 1C) through two strong hydrogen bonds via ASP285 and GLY335 amino acid residues of AChE. Stability of Soy-B and AChE binding was further enhanced by van der Waals' interactions with TYR70, ASP276, VAL277, TRP279, ASN280, SER286, ILE287, TYR334 and LEU358 amino acid residues ( Figure 1D).

Analysis of Interaction of Soy-B with BuChE
Interaction between Soy-B and BuChE revealed that Soy-B binds within the active site cavity of BuChE ( Figure 2). Comparative analysis of Soy-B binding mode was done with native control ligand (Tacrine) and standardization was performed by re-docking the native ligand ( Figure 2C). Hydrophobic pi-pi stacking and pi-alkyl interactions were observed during native ligand binding with BuChE active site, and similar interactions were found when the native ligand was redocked with BuChE ( Figure 2C,E). Not only hydrophobic interactions, but also van der Waals' interactions played an important role in binding of native ligand to the BuChE active site. However, Soy-B showed only van der Waals' interactions via 17 amino acid residues (ASN68, ASP70, GLY78, SER79, TRP82,  GLY116, GLY117, GLN119, THR120, PRO285, SER287, ASN289, ALA328, PHE329,  TYR332, TRP430, HIS438) of BuChE ( Figure 2D).

Analysis of Interaction of Soy-B with BuChE
Interaction between Soy-B and BuChE revealed that Soy-B binds within the active site cavity of BuChE ( Figure 2). Comparative analysis of Soy-B binding mode was done with native control ligand (Tacrine) and standardization was performed by re-docking the native ligand ( Figure 2C). Hydrophobic pi-pi stacking and pi-alkyl interactions were observed during native ligand binding with BuChE active site, and similar interactions were found when the native ligand was redocked with BuChE ( Figure 2C,E). Not only hydrophobic interactions, but also van der Waals' interactions played an important role in binding of native ligand to the BuChE active site. However, Soy-B showed only van der Waals' interactions via 17 amino acid residues (ASN68, ASP70, GLY78, SER79, TRP82,  GLY116, GLY117, GLN119, THR120, PRO285, SER287, ASN289, ALA328, PHE329, TYR332,  TRP430, HIS438) of BuChE ( Figure 2D).

Analysis of Interaction of Soy-B with GS3Kβ
Soy-B interaction with GS3Kβ is shown in Figure 3. The binding mode of Soy-B interaction was compared with native ligand control, ADP ( Figure 3B). In this study, interaction GS3Kβ inhibitor bromoindirubin was also observed. Interestingly, all of the tested compounds, inhibitors and native ligands bound to the same position in the active site cavity of GS3Kβ ( Figure 3B). Superimposition of the native ligand, the redocked native ligand, and the inhibitor at the active site confirmed the standardization of the protocol ( Figure 3C,D). However, hydrogen bonds, electrostatic and van der Waals' interactions dominated the stability of native ligand and GS3Kβ binding ( Figure 3F). On the other hand, Soy-B showed one hydrogen bond with GS3Kβ through GLN185 amino acid residue and one electrostatic interaction with LYS183, followed by seventeen hydrophobic interactions ( Figure 3E).

Analysis of Interaction of Soy-B with GS3Kβ
Soy-B interaction with GS3Kβ is shown in Figure 3. The binding mode of Soy-B interaction was compared with native ligand control, ADP ( Figure 3B). In this study, interaction GS3Kβ inhibitor bromoindirubin was also observed. Interestingly, all of the tested compounds, inhibitors and native ligands bound to the same position in the active site cavity of GS3Kβ ( Figure 3B). Superimposition of the native ligand, the redocked native ligand, and the inhibitor at the active site confirmed the standardization of the protocol ( Figure 3C,D). However, hydrogen bonds, electrostatic and van der Waals' interactions dominated the stability of native ligand and GS3Kβ binding ( Figure 3F). On the other hand, Soy-B showed one hydrogen bond with GS3Kβ through GLN185 amino acid residue and one electrostatic interaction with LYS183, followed by seventeen hydrophobic interactions ( Figure 3E).

Analysis of Interaction of Soy-B with MOA-A
Interaction between Soy-B and MOA-A is shown in Figure 4. Binding mode of Soy-B interaction with MOA-A was compared with the native ligand control, Harmine ( Figure 4A,B). All of the tested compounds and the native ligand bound to the same position in the active site cavity of MOA-A ( Figure 4B). Superimposition of the native ligand and the redocked native ligand at the active site confirmed the standardization of the protocol ( Figure 4C). Redocking of the native ligand ( Figure 4E) showed hydrogen bonds and hydrophobic interactions, whereas Soy-B binding with MOA-A showed two hydrogen bonds with CYS323 and THR336, one electrostatic interaction with PHE352, and 19 hydrophobic interactions ( Figure 4D). cavity of GS3Kβ ( Figure 3B). Superimposition of the native ligand, the redocked native ligand, and the inhibitor at the active site confirmed the standardization of the protocol ( Figure 3C,D). However, hydrogen bonds, electrostatic and van der Waals' interactions dominated the stability of native ligand and GS3Kβ binding ( Figure 3F). On the other hand, Soy-B showed one hydrogen bond with GS3Kβ through GLN185 amino acid residue and one electrostatic interaction with LYS183, followed by seventeen hydrophobic interactions ( Figure 3E).   Figure 4B). Superimposition of the native ligand and the redocked native ligand at the active site confirmed the standardization of the protocol ( Figure 4C). Redocking of the native ligand ( Figure 4E) showed hydrogen bonds and hydrophobic interactions, whereas Soy-B binding with MOA-A showed two hydrogen bonds with CYS323 and THR336, one electrostatic interaction with PHE352, and 19 hydrophobic interactions ( Figure 4D).  Figure 5A,B). The interaction pattern of the best hit (Soyasapogenol B) with the target protein is represented in Figure 5. The protocol was validated through redocking of the native ligand, and it was observed that redocked ligand occupied the same residues (LEU171, CYS172, ILE199, GLN206, ILE316, TYR326) as the native ligand ( Figure 5C,E). It was also observed that Soy-B and MOA-B complexes were stabilized by three (LEU171, ILE198, TYR398) important residues of the catalytic site. Moreover, several van der Waals' interactions were also observed in stabilizing the Soy-B and MOA-B complexes ( Figure 5D).  Figure 5A,B). The interaction pattern of the best hit (Soyasapogenol B) with the target protein is represented in Figure 5. The protocol was validated through redocking of the native ligand, and it was observed that redocked ligand occupied the same residues (LEU171, CYS172, ILE199, GLN206, ILE316, TYR326) as the native ligand ( Figure 5C,E). It was also observed that Soy-B and MOA-B complexes were stabilized by three (LEU171, ILE198, TYR398) important residues of the catalytic site. Moreover, several van der Waals' interactions were also observed in stabilizing the Soy-B and MOA-B complexes ( Figure 5D).

Analysis of Interaction of Soy-B with NMDA
The interaction binding pattern of compounds and native ligands revealed that all of the ligands occupied the active site gorge residues ( Figure 6A-C). Redocked native ligand bound to the same location as of co-crystallized native ligand in NMDA with hydrogen bond (PRO124, THR126, ARG131), electrostatic interaction (ASP224), hydrophobic interactions (PHE16, PHE92, PRO124, VAL227, PHE250), and several van der Waals' interactions ( Figure 6E). The complex of Soy-B and NMDA was stabilized mostly by the van der Waals' interactions and one residue (THR126) were observed to interacted with hydrogen bond ( Figure 6D).

Analysis of Interaction of Soy-B with NMDA
The interaction binding pattern of compounds and native ligands revealed that all of the ligands occupied the active site gorge residues ( Figure 6A-C). Redocked native ligand bound to the same location as of co-crystallized native ligand in NMDA with hydrogen bond (PRO124, THR126, ARG131), electrostatic interaction (ASP224), hydrophobic interactions (PHE16, PHE92, PRO124, VAL227, PHE250), and several van der Waals' interactions ( Figure 6E). The complex of Soy-B and NMDA was stabilized mostly by the van der Waals' interactions and one residue (THR126) were observed to interacted with hydrogen bond ( Figure 6D).

Root Mean Square Deviation (RMSD) Analysis
MD simulation analysis was carried out to understand the selected proteins and Soy-B complexes' nature of interaction (dynamic) and stability. Figure 7 presents the MD simulation of Protein-Soy-B complex for 100 ns. The root-mean-square deviation (RMSD) is a measure that accounts for the protein stability and its deviation from the initial structure during the MD simulation [28]. For the early 10 ns, the AChE, RMSD alone ranged between 0.00 and 2.04 Å. Then, the subsequent simulation was constant. It showed that the AChE structure during the simulation was not considerably altered. For the first 10 ns, and then 0.50 Å for the rest of the simulation, the AChE-Soy-B complex's RMSD varied within 0.00-2.96 Å. The average RMSD values of AChE alone and in complex with Soy-B during 10-100 ns simulation were estimated as 2.18 and 1.94 Å, respectively (Figure 7a). the ligands occupied the active site gorge residues ( Figure 6A-C). Redocked native ligand bound to the same location as of co-crystallized native ligand in NMDA with hydrogen bond (PRO124, THR126, ARG131), electrostatic interaction (ASP224), hydrophobic interactions (PHE16, PHE92, PRO124, VAL227, PHE250), and several van der Waals' interactions ( Figure 6E). The complex of Soy-B and NMDA was stabilized mostly by the van der Waals' interactions and one residue (THR126) were observed to interacted with hydrogen bond ( Figure 6D).  MD simulation analysis was carried out to understand the selected proteins and Soy-B complexes' nature of interaction (dynamic) and stability. Figure 7 presents the MD simulation of Protein-Soy-B complex for 100 ns. The root-mean-square deviation (RMSD) is a measure that accounts for the protein stability and its deviation from the initial structure during the MD simulation [28]. For the early 10 ns, the AChE, RMSD alone ranged between 0.00 and 2.04 Å. Then, the subsequent simulation was constant. It showed that the AChE structure during the simulation was not considerably altered. For the first 10 ns, and then 0.50 Å for the rest of the simulation, the AChE-Soy-B complex's RMSD varied within 0.00-2.96 Å. The average RMSD values of AChE alone and in complex with Soy-B during 10-100 ns simulation were estimated as 2.18 and 1.94 Å, respectively (Figure 7a).

Protein-Ligand Interaction Analysis
As seen in Figure 11, hydrogen bonds and hydrophobic interactions play an essential role in sustaining and stabilizing the complexes of targeted proteins and Soy-B. In addition, water bridges and ionic association are also critical to forming a stable protein-ligand complex. The amino acid residues SER235, GLU306, and ASN525, were the most important in forming hydrogen bonds between AChE and Soy-B (Figure 11a). Similarly, amino acid residues such as ASN68, TRP82, and ASN289 were observed to stabilize the BuChE-Soy-B complex through hydrogen bonding, while PRO285, ALA328, PHE329 and TYR332 were involved in hydrophobic interactions (Figure 11b). Likewise, in the formation of the MAO-A and Soy-B complex, amino acid residues such as ASN125, ASN212, and GLU485 were involved in hydrogen bond formation, while ALA111, TYR124, and TRP128 were engaged with Soy-B through hydrophobic interactions (Figure 11c). Moreover, in MAO-B and Soy-B complex formation, amino acid residue ARG307 was predominantly involved in hydrogen bond formation, and TYR80 formed hydrophobic interaction (Figure 11d). In addition, GS3Kβ and Soy-B complex formed amino acid residues VAL135 and LYS183 hydrogen bond formation while ALA83, LEU132, and LEU188 formed hydrophobic interactions (Figure 11e). Moreover, in NMDA and Soy-B complex formation, amino acid residue SER248, PHE246, ASP224, and GLN144 were predominantly involved in hydrogen bond formation, and TYR184, ILE183, and LEU146 formed hydrophobic interactions ( Figure 11f).

Protein-Ligand Interaction Analysis
As seen in Figure 11, hydrogen bonds and hydrophobic interactions play an essential role in sustaining and stabilizing the complexes of targeted proteins and Soy-B. In addition, water bridges and ionic association are also critical to forming a stable protein-ligand complex. The amino acid residues SER235, GLU306, and ASN525, were the most important in forming hydrogen bonds between AChE and Soy-B (Figure 11a). Similarly, amino acid residues such as ASN68, TRP82, and ASN289 were observed to stabilize the BuChE-Soy-B complex through hydrogen bonding, while PRO285, ALA328, PHE329 and TYR332 were involved in hydrophobic interactions (Figure 11b). Likewise, in the formation of the MAO-A and Soy-B complex, amino acid residues such as ASN125, ASN212, and GLU485 were involved in hydrogen bond formation, while ALA111, TYR124, and TRP128 were engaged with Soy-B through hydrophobic interactions (Figure 11c). Moreover, in MAO-B and Soy-B complex formation, amino acid residue ARG307 was predominantly involved in hydrogen bond formation, and TYR80 formed hydrophobic interaction (Figure 11d). In addition, GS3Kβ and Soy-B complex formed amino acid residues VAL135 and LYS183 hydrogen bond formation while ALA83, LEU132, and LEU188 formed hydrophobic interactions (Figure 11e). Moreover, in NMDA and Soy-B complex formation, amino acid residue SER248, PHE246, ASP224, and GLN144 were predominantly involved in hydrogen bond formation, and TYR184, ILE183, and LEU146 formed hydrophobic interactions (Figure 11f).

Protein-Ligand Interaction Analysis
As seen in Figure 11, hydrogen bonds and hydrophobic interactions play an essential role in sustaining and stabilizing the complexes of targeted proteins and Soy-B. In addition, water bridges and ionic association are also critical to forming a stable protein-ligand complex. The amino acid residues SER235, GLU306, and ASN525, were the most important in forming hydrogen bonds between AChE and Soy-B (Figure 11a). Similarly, amino acid residues such as ASN68, TRP82, and ASN289 were observed to stabilize the BuChE-Soy-B complex through hydrogen bonding, while PRO285, ALA328, PHE329 and TYR332 were involved in hydrophobic interactions (Figure 11b). Likewise, in the formation of the MAO-A and Soy-B complex, amino acid residues such as ASN125, ASN212, and GLU485 were involved in hydrogen bond formation, while ALA111, TYR124, and TRP128 were engaged with Soy-B through hydrophobic interactions (Figure 11c). Moreover, in MAO-B and Soy-B complex formation, amino acid residue ARG307 was predominantly involved in hydrogen bond formation, and TYR80 formed hydrophobic interaction (Figure 11d). In addition, GS3Kβ and Soy-B complex formed amino acid residues VAL135 and LYS183 hydrogen bond formation while ALA83, LEU132, and LEU188 formed hydrophobic interactions (Figure 11e). Moreover, in NMDA and Soy-B complex formation, amino acid residue SER248, PHE246, ASP224, and GLN144 were predominantly involved in hydrogen bond formation, and TYR184, ILE183, and LEU146 formed hydrophobic interactions (Figure 11f). The involvement of amino acid residues in terms of percent simulation time during the formation of a complex between Soy-B and targeted proteins (AChE, BuChE, MAO-A, MAO-B, GS3Kβ and NMDA) was also established ( Figure 12). Only those amino acid residues which interacted with ligand for more than 30% simulation time are shown. It is clear that during the simulation of ACHE-Soy-B, the amino acid residues SER235, LEU305, and ASN525 were involved with the ligand for about 36%, 43% and 41% of simulation time, respectively (Figure 12a). Similarly, during the BuChE-Soy-B simulation, TRP82 and ASN289 were engaged for 96% and 68%, respectively (Figure 12b). Likewise, ASN125 of MAO-A interacted with Soy-B for about 45% of simulation time through a water bridge (Figure 12c), while ARG307 of MAO-B formed two interactions with two different hydroxyl groups of Soy-B for 59% and 35%, respectively (Figure 12d). Soy-B formed interactions with GS3Kβ at the VAL135 amino residue, forming hydroxyl bonds for about 76% of the simulation (Figure 12e), whereas with NMDA, hydroxyl interactions were formed with the target proteins PHE246, SER248, and ASP224, which were stable throughout 94%, 87% and 40% of the simulation time, respectively. The involvement of amino acid residues in terms of percent simulation time dur the formation of a complex between Soy-B and targeted proteins (AChE, BuChE, MA A, MAO-B, GS3Kβ and NMDA) was also established ( Figure 12). Only those amino a residues which interacted with ligand for more than 30% simulation time are shown. I clear that during the simulation of ACHE-Soy-B, the amino acid residues SER235, LEU3 and ASN525 were involved with the ligand for about 36%, 43% and 41% of simulat time, respectively (Figure 12a). Similarly, during the BuChE-Soy-B simulation, TRP82 a ASN289 were engaged for 96% and 68%, respectively (Figure 12b). Likewise, ASN125 MAO-A interacted with Soy-B for about 45% of simulation time through a water brid (Figure 12c), while ARG307 of MAO-B formed two interactions with two different h droxyl groups of Soy-B for 59% and 35%, respectively (Figure 12d). Soy-B formed inter tions with GS3Kβ at the VAL135 amino residue, forming hydroxyl bonds for about 7 of the simulation (Figure 12e), whereas with NMDA, hydroxyl interactions were form with the target proteins PHE246, SER248, and ASP224, which were stable throughout 94 87% and 40% of the simulation time, respectively. In addition, the formation of a stable protein and ligand complex was established determining the total number of contacts formed between them during simulation (Figu 13). It is clear that during the simulation, the total number of contacts between Soy-B a AChE, BuChE, MAO-A, MAO-B, GS3Kβ, and NMDA varied between 0 and 8, 2 and 9 and 9, 0 and 7, 2 and 9, and 3 and 9, respectively. On average, AChE, BuChE, MAO MAO-B, GS3Kβ and NMDA formed 5, 5, 4, 4, 5 and 5 contacts with Soy-B. In addition, the formation of a stable protein and ligand complex was established by determining the total number of contacts formed between them during simulation ( Figure 13). It is clear that during the simulation, the total number of contacts between Soy-B and AChE, BuChE, MAO-A, MAO-B, GS3Kβ, and NMDA varied between 0 and 8, 2 and 9, 0 and 9, 0 and 7, 2 and 9, and 3 and 9, respectively. On average, AChE, BuChE, MAO-A, MAO-B, GS3Kβ and NMDA formed 5, 5, 4, 4, 5 and 5 contacts with Soy-B. Thus, from the results it can be inferred that the Soy-B was shown to have good affinity, and tends to form stable complexes with hydroxyl group interactions.

Discussion
Alzheimer's Disease (AD) is a neurodegenerative disease affecting 55 million people globally, which is expected to increase over 106.4 million by 2050 [5,29]. Even though there are tremendous research efforts directed toward AD treatment, an effective drug is still lacking. Therefore, there is a need for alternate strategies for AD treatment. Since phytochemicals remain a partially explored pool of therapeutic strategies, the present study was intended to virtually screen the 44 natural BDNF inducers against the AChE, BuChE, GSK3β, MOAA, MOAB, NDMA proteins to identify potential therapeutic candidates for AD treatment. AChE is a key enzyme of the cholinergic system, and has been reported to be deteriorated in the progression of AD, which is accompanied by the decline of acetylcholine and the loss of cholinergic neurons in the forebrain. There are multiple therapeutic modalities targeting the cholinergic system; in particular, AChE is the main target that has been proven to be a promising therapeutic strategy [30]. BuChE is a serine hydrolase enzyme that plays a crucial role in the regulation of acetylcholine metabolism, and is closely related to AChE. It tends to be predominantly expressed in the white matter, glia and neurons regulating cognition and behavior in the AD condition [31]. In an AD mouse model, accumulation of fibrillary Aβ plaques indicated an active role in the pathological development of AD. Using a strong and specific BuChE inhibitor to increase acetylcholine levels and intimidate fibrillar Aβ accumulation has been demonstrated to be an easy and successful treatment strategy for AD [32].
Both Monoamine oxidase-A (MAO-A) and MAO-B are involved in AD pathogenesis. Neurodegeneration in AD is related to the altered function of any isoform or any associated disturbance of neurotransmitter substrates, like serotonin, dopamine, or noradrenaline. MAO might possibly play a role in neuropathology due to the production of hydrogen peroxide as a byproduct of the deamination reaction. When antioxidant mechanisms are weakened, such as during aging, and especially in AD, the resulting oxidative stress could lead to cell death, which inevitably involves the role of mitochondria. MAO-A-immunoreactive cell loss is amplified in late-stage cognitive impairment in brainstem monoaminergic nuclei and other locations. MAO-B is predominantly expressed in glia and has been thoroughly investigated with respect to its function in neurodegeneration. MAO-A/B-associated shift in AD levels of aminergic neurotransmitters likely contributes to the Thus, from the results it can be inferred that the Soy-B was shown to have good affinity, and tends to form stable complexes with hydroxyl group interactions.

Discussion
Alzheimer's Disease (AD) is a neurodegenerative disease affecting 55 million people globally, which is expected to increase over 106.4 million by 2050 [5,29]. Even though there are tremendous research efforts directed toward AD treatment, an effective drug is still lacking. Therefore, there is a need for alternate strategies for AD treatment. Since phytochemicals remain a partially explored pool of therapeutic strategies, the present study was intended to virtually screen the 44 natural BDNF inducers against the AChE, BuChE, GSK3β, MOAA, MOAB, NDMA proteins to identify potential therapeutic candidates for AD treatment. AChE is a key enzyme of the cholinergic system, and has been reported to be deteriorated in the progression of AD, which is accompanied by the decline of acetylcholine and the loss of cholinergic neurons in the forebrain. There are multiple therapeutic modalities targeting the cholinergic system; in particular, AChE is the main target that has been proven to be a promising therapeutic strategy [30]. BuChE is a serine hydrolase enzyme that plays a crucial role in the regulation of acetylcholine metabolism, and is closely related to AChE. It tends to be predominantly expressed in the white matter, glia and neurons regulating cognition and behavior in the AD condition [31]. In an AD mouse model, accumulation of fibrillary Aβ plaques indicated an active role in the pathological development of AD. Using a strong and specific BuChE inhibitor to increase acetylcholine levels and intimidate fibrillar Aβ accumulation has been demonstrated to be an easy and successful treatment strategy for AD [32].
Both Monoamine oxidase-A (MAO-A) and MAO-B are involved in AD pathogenesis. Neurodegeneration in AD is related to the altered function of any isoform or any associated disturbance of neurotransmitter substrates, like serotonin, dopamine, or noradrenaline. MAO might possibly play a role in neuropathology due to the production of hydrogen peroxide as a byproduct of the deamination reaction. When antioxidant mechanisms are weakened, such as during aging, and especially in AD, the resulting oxidative stress could lead to cell death, which inevitably involves the role of mitochondria. MAO-Aimmunoreactive cell loss is amplified in late-stage cognitive impairment in brainstem monoaminergic nuclei and other locations. MAO-B is predominantly expressed in glia and has been thoroughly investigated with respect to its function in neurodegeneration. MAO-A/B-associated shift in AD levels of aminergic neurotransmitters likely contributes to the neurobiology of a spectrum of neurological disorders in AD populations [33]. GSK3β leads to tau protein hyperphosphorylation, the major component of neurofibrillary tangles (NFTs), which are hallmarks of AD. In addition, GSK3β regulates several neuronal functions, such as Aβ-induced cell death, axonal transportation, cholinergic functioning, and adult neurogenesis that are dysregularized throughout the pathogenesis of AD [34]. NMDA plays a key role in the synaptic transmission and plasticity of learning, memory, and the development of the nervous system, as well as in neurotoxicity. Ca2+ signaling causes neuronal cell death, which is clinically linked to cognitive impairment and the establishment of pathological architecture in AD. This also justifies the clinical trial of the NMDA receptor antipathist memantine as a protective and symptomatologic treatment for AD [35].
A total of 44 natural BDNF inducer compounds were assessed for their toxicity, and the physicochemical properties of all of these compounds were evaluated. TPSA was analyzed in order to assess the efficiency of crossing BBB, and it was found that Imperatorin, Eugenol, Auraptene, Soyasapogenol B (Soy-B), Tanshinone II-A, Caffeine, Honokiol and Magnolol possessed TPSA < 60. The BBB is a characteristic barrier that limits the efficiency of therapeutic modalities by reducing their bioavailability as a result of their inability to cross the BBB [36]. Therefore, crossing the BBB is crucial in treating neurodegenerative disorders. Cumulative analysis of the Lipinski rule and toxicity in the BDNF inducers showed that Soy-B, Huperzine A, Hydroxytyrosol, Alpinetin, and Scopoletin demonstrated good absorption efficiency without violating the Lipinski rule or exhibiting toxicity.
Molecular docking analysis suggested that Soy-B had a binding energy score of about −8.7 Kcal/mol, −9.9 Kcal/mol, −8.7 Kcal/mol, −7.9 Kcal/mol, −7.6 Kcal/mol and −8.4 Kcal/mol against the target proteins AChE, BuChE, GSK3β, MOAA, MOAB and NDMA, respectively, on the basis of which Soy-B was selected as a candidate for further simulation studies.
Soy-B interacted with AChE through TYR70, ASP276, VAL277, TRP279, ASN280, ASP285, SER286, ILE287, TYR334, GLY335 and LEU358 amino acid residues ( Figure 1D). In another study, semi-synthetic piperidine alkaloids also interacted with the same amino acid residues TYR70, ASP285, SER286, TYR334 and GLY335 of the active site of AChE [37]. On the other hand, the interaction of Soy-B with BuChE showed the involvement of ASN68, ASP70, GLY78, SER79, TRP82, GLY116, GLY117, GLN119, THR120, PRO285, SER287, ASN289, ALA328, PHE329, TYR332, TRP430, and HIS438 amino acid residues ( Figure 2D). Similary, all of these amino acids were reported to be involved in the interaction of Bacoside X with BuChE in a recent study [38]. Soy-B exhibited an interaction with the ATP binding pocket of GS3Kβ through GLN185 and LYS183 amino acid residues ( Figure 3E). In a similar way, Salvianolic acid B interacted with the same active site pocket of GS3Kβ via the GLN185 and LYS183 amino acids [39]. Moreover, Soy-B binding with MOA-A showed a strong interaction via the CYS323, THR336, and PHE352 amino acid residues ( Figure 4D), while Soy-B and MOA-B interaction was stabilized by the LEU171, ILE198, and TYR398 amino acid residues ( Figure 5D). Thiol of CYS323 of MOA-A plays a crucial role in interaction with different compounds, such as piperine and β -carboline analogs [40,41]. Importantly, Soy-B also showed interaction with the CYS323 residue of MOA-A. In one study, desmodeleganine (a natural alkaloid) was shown to interact with the LEU171 and ILE198 amino acid residues of MOA-B [42]. Meanwhile, several natural compound derivatives interacted with the aromatic cage of MOA-B via TYR398 [43]. In addition, inhibitor binding and the important role of the TYR398 amino acid residue of MOA-B were also reported in a recent study [44]. Interestingly, Soy-B was shown to interact with all of these crucial amino acid residues of MOA-B. Soy-B and NMDA interaction was stabilized by strong bonding with THR126 ( Figure 6D). Similarly, catechin, gingetin and ginkgolic acid showed strong hydrogen bonding with THR126 [45]. Furthermore, MS analysis (RMSD analysis, RMSF analysis, secondary structure, and protein ligand interaction analysis) during 0-100 ns indicated that the Soy-B had formed stable complexes with the targeted proteins (AChE, BuChE MAO-A, MAO-B, GS3Kβ, and NMDA). Figure 7a-f presents the average RMSD value for Soy-B -protein complex for AChE, BuChE, GSK3β, MOAA, MOAB and NDMA as 1.94 Å, 2.11 Å, 5.07 Å, 2.56 Å, 3.83 Å and 4.07 Å, respectively. However, none of the RMSD values fluctuated more than the acceptable limit of 2 Å [46].
Soy-B belongs to the soyasaponins predominantly isolated from Glycine max and shown to exert various pharmacological properties (anti-inflammatory [47], phytooestrogenic [48] and memory enhancement [49,50]). Soy-B has been shown to attenuate LPSinduced memory impairment in mice via NF-κB-mediated BDNF expression and prevent scopolamine-induced memory impairment [50,51]. In LPS-treated mice and corticosteronestimulated SH-SY5Y cells, Soy-B enhanced cAMP response element-binding protein phosphorylation and BDNF expression, and reduced NF-B activation. Soyasaponin subtypes A1, A2 and I reduced LPS-induced cyclooxygenase 2 (COX-2) production via inhibiting NF-κB. Soyasaponins appear to reduce memory loss by inhibiting NF-B-mediated inflammation, according to this research [52,53]. In the present study, Soy-B was predicted to be the best among all of the compounds screened. However, the limitation of the present study is that our findings are only based on the in silico approach, and need to be further confirmed via in vitro and in vivo experiments.

Calculation of Physicochemical Properties and Prediction of Toxicity Potential
A total of 44 Natural BDNF inducer compounds were assessed for their toxicity, and the physicochemical properties of all of the compounds were assessed using the Orisis Datawarrior property explorer tool [54]. The tool was also used to assess the molecular weight, the number of hydrogen bond donors and acceptors, rotatable bonds, cLogP value, topological polar surface area, and Lipinski's rule violation [55]. The absorption % was calculated per Zhao et al. [56] using the formula Absorption % = 109 − (0.345 × TPSA). In the Orisis Datawarrior tool's analysis, the predicted toxicity values were dependent on comparing the precalculated investigated molecules with the structures of the tested molecules. Various aspects and effects of the toxicity, including tumorigenicity, mutagenicity, and irritability of the tested compounds, were also tested using the Orisis Datawarrior tool.  (PubChem ID:131682), respectively, were retrieved from the Protein Data Bank database. Each target was prepared using the AutoDock 4.2 tool, which included solvation parameters, polar hydrogen and Kollman united atom charges and converted the target into pdbqt format.

Ligand Preparation
The three-dimensional structure of ligands was retrieved from the PubChem database. OpenBabel tool (http://openbabel.org (accessed on 21 December 2021)) was used to convert each ligand into pdbqt format [57].

Molecular Docking
Docking studies were carried out using AutoDock Vina [58]. For each target on the specified active site, a grid box was established before docking. The grid configurations of x, y, and z coordinates were set to 60 × 60 × 60, with the grid box's center set to x: 20.33, y: 16.53 z: −10.38 for GSK3β. The center of x: 4.468, y: 70.009, z: 65.930 for AChE was selected through discovery studio visualizer (BIOVIA) from the attributes of docked native ligand in its specific target protein, while the center was set as x: 5.61, y: 38, z: −17.16 for NMDA [14,59]. The grid configurations of x, y, and z coordinates were set to 33 × 33 × 33 with the grid box center set to x: 140.117, y: 122.247, z: 38.986 for BuChE [60]. Meanwhile, the grid configurations of x, y, and z coordinates were set to 30 × 30 × 30 with grid box center set as x: 40.54, y: 27.04, z: −14.49 for MAO-A and x: 52.34, y: 156.14, z: 28.34 for MAO-B as discussed in previous reports [14,61]. The AutoDock Vina results are displayed as an affinity (kcal/mol). For each ligand-target interaction, the algorithm is divided into 10 modes in descending order. With the PyMol tool's aid, the best conformation with the lowest energy was chosen and saved as a complex. However, using Discovery Studio Visualizer, the complex produced was further visualized and analyzed (BIOVIA, San Diego, CA, USA, Dassault Systèmes, Velizy, France) [60]. The binding affinity (Kd) of ligands for the target protein was calculated from the Gibb's free binding energy (∆G) using the following relation: In this equation, the universal gas constant is denoted as "R", whereas temperature is defined as "T".
The ligands with the minimum Gibb's free binding energy were selected for further analysis.

Molecular Dynamic (MD) Simulation Study
The stability and dynamic nature of the ligand inside a protein's binding pocked were determined by performing MD simulation for 100 ns using Desmond (Schrödinger, LLC, New York, NY, USA), as described previously [62]. The initial structures of the proteinligand complex for MD simulation were obtained from the docking pose. Protein-ligand complex was prepared using Protein Preparation wizard in the Maestro module, which included complex optimization and minimization. The system creator tool prepared all structures within the orthorhombic simulation box, which was then solvated with TIP3P explicit water molecules. The models were neutralized by inserting appropriate counter ions and 0.15 M salt (NaCl) for the physiological condition simulation. LSSB-2005 was used to balance the system with 2000 iterating steps, and forcefield was used to maintain the convergence 1 kcal/mol/Å criteria. MD simulation was performed for 100 ns with an NPT set at 300 K temperature with 1 atm pressure. A Nose-Hoover thermostat and Matrtyna-Tobias-Klein barostat were used to maintain these parameters [63,64]. During simulation, a time step of 2 fs was established, and the energy/structures were recorded every 10 ps. To establish protein-lining complex stability, the metrics like Root-Medium-Square Deviation (RMSD), Root-Medium-Square-Fluctuation (RMSF), radius of gyration (Rg) and Solvent-Accessible Surface Area (SASA) were examined.

Conclusions
The present study conducted a screening of potent abrineurin/BDNF natural inducers to find multi-targeting agents against Alzheimer's Disease. After rigorous initial in silico screening via physicochemical properties, the Lipinski rule, total polar surface area, toxicity assessment, and oral absorption efficiency, five compounds (Soy-B, Huperzine A, Hydroxytyrosol, Alpinetin and Calycosin) were selected. Further molecular docking and interaction analysis against different Alzheimer's targets (AChE, BuChE MAO-A, MAO-B, GS3Kβ and NMDA) predicted Soy-B to be the most promising compound. However, the stability of Soy-B binding with the target proteins was established by molecular dynamic simulation. Simulation analysis revealed that Soy-B interaction with the active site amino acid residues of the target proteins retained for most of the simulation time, which forms a stable complex. Further in vitro and in vivo explorations are needed to establish Soy-B as a therapeutic candidate for Alzheimer's disease. However, it is noteworthy to mention that BDNF induction could effectively mitigate neurological ailments, and exploring a solution on BDNF inducer scaffolds for Alzheimer's appears reasonable. Nevertheless, the present preliminary study predicted that the Soy-B scaffold has all of the characteristics needed for multi-targeting solutions against complex disease like Alzheimer's.